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1 Introduction 

Random matrices play an important role in physics and mathematics ISni US El 13123 El US- it 


has been observed more and more in the recent years how deeply random matrices are related to 


integrability (r-functions), and algebraic geometry. 

Here, we consider the computation of large n asymptotics for orhogonal polynomials as an example 
of a problem where the concepts of integrability, isomonodromy and algebraic geometry appear and 
combine. 

The method presented here below, is not, to that date, rigorous mathematicaly. It is based on the 
asumption that an integral with a large number of variables can be approximated by a saddle-point 
method. This asumption was never proven rigorously, it is mostly based on “physical intuition”. 
However, the results given by that method have been rigorously proven by another method, namely 
the Riemann-Hilbert method 0013112] The method presented below was presented in many 


works [3130 EDI d. 


2 Definitions 


Here we consider the 1-Hermitean matrix model with polynomial potential: 



( 2 - 1 ) 


^E-mail: 



where A(xi,... ,xn) ■= ~ ^j)y eigenvalues of the matrix M, and V{x) 

is a polynomial called the potential: 

deg V 

V{x) = ^ gkx’^ (2-2) 

k=0 

Remark 2.1 All the calculations which are presented below, can be extended to a more general 
setting, with no big fundamental changes: 

- one can consider V\x) any rational fraction H instead of polynomial, in particular one can 
add logarithmic terms to the potential V(x). 

- one can consider arbitrary paths (or homology class of paths) of integrations insteaf ofM.^, 
in particular finite segments H 

- one can study non hermitean matrix models m, where the Vandermonde is replaced by A^ 
where P = 1,2, 4. 

- one can consider multi-matrix models, in particular 2-matrix model 


3 Orthogonal polynomials 


Consider the family of monic polynomials Pn{x) = x'^+0{x'^ dehned by the orthogonality relation: 




(3-3) 


It is well known that the partition function is given by El: 


7V-1 

Zn = N\1[ hn 

n=0 


(3-4) 


Such an orthogonal family always exists if the integration path is M or a subset of M, and if 
the potential is a real polynomial. In the more general setting, the orthogonal polynomials “nearly 
always” exist (for arbitrary potentials, the set of paths for which they don’t exist is enumerable). 


We dehne the kernel: 

K{x,y) 

One has the following usefull theorems: 


N-l 

E 

n=0 


Pn{x)pn{y) 

hri. 


(3-5) 


Theorem 3.1 Dyson’s theorem m-- any correlation function of eigenvalues, can be written in terms 
of the kernel K: 


p(A,,.... At) = det(A'(Aj, A,)) 


(3-6) 



Thus, if one knows the orthogonal polynomials, then one knows all the correlation functions. 


Theorem 3.2 Christoffel-Darboux theorem 123 Ef.- The kernel K{x,y) can be written: 

PN{.x)pN-i{y) - PN{.y)PN-i{x) 


K{x,y) = 77V 


x-y 


(3-7) 


Thus, if one knows the polynomials pn and Pn- 1 i then one knows all the correlation functions. 

Our goal now, is to hnd large N ’’strong” asymptotics for pjv and PN-h in order to have the large 
N behaviours of any correlation functions. 

Notation: we dehne the wave functions: 

V’n(a:) := (3-8) 

Vhn 

they are orthonormal: 

I ^nm (^”9) 


4 Differential equations and integrability 


It can be proven that {'ipn,4’n-i) obey a differential equation of the form [313 EHl EHl E] : 


J_ y / 'lpn{x) 
N dx V'0n-i(a:) 


Vn{x) 


f i^n{x) \ 

\^/Jn-l{x) ) 


(4-10) 


where T>n{x) is a 2 x 2 matrix, whose coefficients are polynomial in x, of degree at most deg Id'. (In 
case V is a rational function, then P is a rational function with the same poles). 

(V’n, V’n-i) also obeys differential equations with respect to the parameters of the model [313, he. 
the coupling constants, i.e. the gus defined in 12-21 


i^n{x) \ 

N dgk y'llJn-ii.x)) 



(4-11) 


where Un,k{x) is a 2 x 2 matrix, whose coefficients are polynomial in x, of degree at most k. 

It is also possible to find some discrete recursion relation in n (see jS])- 

The compatibility of these differential systems, i.e. as well as 

J -J ■> dx dgk ogk ax ’ dgj dgk agk agj ’ 

compatibility with the discrete recursion, imply integrability, and allows to define a r-function 


123 E]- 


We define the spectral curve as the locus of eigenvalues of Vn{x): 


En{x, y) := det(|/l - Vn{x)) 


(4-12) 





Remark 4.1 In the 1-hermitean-matrix model, is a 2 x 2 matrix, and thus degy En{x,y) = 2, i.e. the 
curve En{x,y) = 0 is an hyperelliptical curve. In other matrix models, one gets algebraic curves which 
are not hyperelliptical. 

Remark 4.2 What we will se below, is that the curve E]\f{x,y) has a large N limit E{x,y), which is also 
an hyperelliptical curve. In general, the matrix T)]\f{x) has no large N limit. 


5 Riemann-Hilbert problems and isomonodromies 


The 2x2 system Vjq has 2 independent solutions: 


'4^n{x) \ ( i’nix) 

N dx \i>n-i{x) ) I r/>n_i(a;) 


N dx I 0n-i(a:) 


= Tdnix) 


where the wronskian is non-vanishing: det ( \ ] 7^ 0- 

We define the matrix of fundamental solutions: 

'tpn{x) (t>n(x) 


^ ■ ' ^n-l(x) (/)n-l(x) 


it obeys the same differential equation: 


lA 

N dx 


tjx) = 'D„(x) fn(x) 


<i>n-l(x) 


(5-13) 


(5-14) 


(5-15) 


Here, the second solution can be constructed explicitely: 


(x) = 


dx' 


X — X‘ 


-'ilJn{x')e 2 




(5-16) 


Notice that (j)n{x) is discontinuous along the integration path of x' (i.e. the real axis in the most 
simple case), the discontinuity is simply 2m'ipn{x). In terms of fundamental solutions, one has the 
jump relation: 

T„(a;-l-z0) = T„(a;-20) (5-17) 

Finding an invertible piecewise analytical matrix, with given large x behaviours, with given jumps 
on the borders between analytical domains, is called a Riemann-Hilbert problem [7||H1I1]. 

It is known that the Riemann-Hilbert problem has a unique solution, and that if two R-H problems 
differ by e (i.e. the difference between jumps and behaviours at oo is bounded by e), then the two 
solutions differ by at most e (roughly speeking, harmonic functions have their extremum on the 
boundaries). Thus, this approach can be used j71 HH |T2] in order to find large N asymptotics of 
orthogonal polynomials: The authors of [Zj considered a guess for the asymptotics, which satisfies 
another R-H problem, which differs from this one by 0(l/iV). 



Notice that the jump matrix in 15-171 is independent of x, of n and of the potential, it is a constant. 
The jump matrix is also called a monodromy, and the fact that the monodromy is a constant, is 
called isomonodromy property EH- 

Consider an invertible, piecewise analytical matrix T„(x), with appropriate behaviours^ at oo, 
which satishes [5-171 then, it is clear that the matrix —has no discontinuity, and 
given its behaviour at oo, it must be a polynomial. Thus, we can prove that T„(x) must satisfy 
a differential system 'Dn{x) with polynomial coefficients. Similarly, the fact that the monodromy is 
independent of Qk and n implies the deformation equations, as well as the discrete recursion relations. 

Thus, the isomonodromy property, implies the existence of compatible differential systems, and 
integrability [nil211Enil2IlE3E]- 

6 WKB—like asymptotics and spectral curve 

Let us look for a formal solution of the form: 

T^(x) = An{x) (6-18) 

where T(x) = diag(Ti(x), T 2 (x)) is a diagonal matrix, and is independent of x. The differential 
system Vn^x) is such that: 

T>Nix) = = An{x)T'{x)A^\x) - ^A'^{x)A^\x) 

= AN{x)T'{x)Ajf\x)+ 0{^) (6-19) 

this means, that, under the asumption that A]\f{x) has a large N limit A{x), T[{x) and TJ^^iyx) are the 
large N limits of the eigenvalues of T>]sf{x). 

With such an hypothesis, one gets for the orthogonal polynomials: 

Tp^ix) ~ ( 6 - 20 ) 

We are now going to show how to derive such a formula. 


7 Orthogonal polynomials as matrix integrals 

7.1 Heine’s formula 


Theorem 7.1 Heine’s theorem 1,V^ . The orthogonal polynomials p — n{x) are given by: 

J dxi... dxN n^i(^ - ^i) (^(a^n • • •, XN)f rij^i 


Pn(0 = 


/ dxi .. . dxN (A(xi,..., xnW n*=i 


^The behaviours at oo are far beyond the scope of this short lecture. They are easily obtained by computing 4>n{x) 
by saddle point method at large x. 








(det(fl - M)) 


(7-21) 


i.e. the orthogonal polynomial is the average of the characteristic polynomial of the random matrix. 

Thus, we can dehne the orthogonal polynomials as matrix integrals, similar to the partition 
function Z dehne in o 


7.2 Another matrix model 

Dehne the potential: 


Vh{x) := V{x) — h In — x) 


and the partition function: 


Zn{h,T):=e t- 




/ » £/ 
dxi... dxn (A(a:i,..., Xn)Y e 

i=i 


-^Vhixi) 


(7-22) 


(7-23) 


i.e. .^Ar(0, 1) = Z is our initial partition function. 
Heine’s formula reads: 


PniO 


^ )-Fn(0,f)) 

^n(0, ^) 


(7-24) 


The idea, is to perform a Taylor expansion in h close to 0 and T close to 1. 


7.2.1 Taylor expansion 


We are interested in n = iV and n = iV — 1, thus T = ^ = 1 + = 1 -|- 0(1/N) and h = 0 or 

h = l/N, i.e. h = 0{l/N): 

T =1 + 0{1/N) , h = 0{l/N) (7-25) 

Roughly speaking: 

~ + 0(1/Ar)) (7-26) 


where all the derivatives are computed at T = 1 and h = 0. 


7.2.2 Topological expansion 

Imagine that Fn has a 1/u^ expansion of the form: 

F = + \f^^^ + O(^) 


(7-27) 






where all and are analytical fnnctions of T and h, than one needs only in order to 
compnte the asymptotics 17-2(11 

Actnaly, that hypothesis is not always trne. It is wrong in the so called ’’mntlicut” case. Bnt it 
can be adapted in that case, we will come back to it in section I11.2I For the moment, let ns condnct 
the calculation only with F^^\ 


8 Computation of derivatives of 


We have dehned: 


this implies that: 


i.e. 


Zn{h,T) ^ 


v? dFr, / n 


T2 dh 


= ^^tr In 


Vh 


(8-28) 

(8-29) 


dFn 

dh 


(8 - 30) 


n 


(tr In {i-M))y^ 


It is a primitive of —^ (tr In [x — M))y^, which behaves as —^ Inx -|- 0(l/x) at large x. Therefore, 
we define the resolvent W{x)\ 


T 1 

W{x) := — ( tr ■ 


n \ X — M 


(8-31) 


Vh 


Notice that it depends on ^ through the potential 14, i.e. through the average < . >. And we dehne 
the effective potential: 


T. 


V^q{x) = Vh{x) - 2T In a; - 2 / {W{x) - -)dx' 


X' 


(8-32) 


which is a primitive of Vl^{x) — 2fF(a;). Thus , we have: 

dFr, 1 


dh 2 


= 7T (r.»(0 - n(0) 


(8-33) 


We also introduce: 


n(x) := , InA(x) := Inx -h f {Q{x') - ^)dx' = ~-^Ves{x) (8-34) 


dT 


H{x,^):= , lnF(0:=/ H{x',^)dx' 


(8-35) 










i.e. 


dh^ 


InHiO 


With these notations, the asymptotics are: 


dhdT 


lnA(0 


MO ~ e-T^eff(C) (i + 0(i/iv)) 


(8-36) 


(8-37) 


Now, we are going to compute W, A, H, etc, in terms of geometric properties of an hyperelliptical 
curve. 


Remark 8.1 This is so far only a sketch of the derivation, valid only in the 1-cut case. In general, has 
no l/v? expansion, and that case will be addressed in section ^^21 

Remark 8.2 These asymptoics are of the form of 16-181 in sectional and thus, ^V'{x) — W{x) is the limit 
of the eigenvalues of ^^^{x). 


9 Saddle point method 

There exists many ways of computing the resolvent and its derivatives with respect to h, T, or other 
parameters. The loop equation method is a very good method, but there is not enough time to 
present it here. There are several saddle-point methods, which all coincide to leading order. We 
are going to present one of them, very intuitive, but not very rigorous on a mathematical ground, 
and not very appropriate for next to leading computations. However, it gives the correct answer to 
leading order. 


Write: 


where the action is: 


ZniK T) = ^ 


5(xi,..., x„) := — 'V Vh{xi) - 2— In (xi - Xj) 

2=1 i>j 


(9-38) 


(9-39) 


The saddle point method consists in hnding conhgurations x* = Xi where S is extremal, i.e. 


Vi = 1,... n, 

i.e., we have the saddle point equation: 


dxi 


= 0 


(9-40) 


tv 7 


Vi = 1,... n. 


(9-41) 










The saddle point approximation^ consists in writting: 

1 


Zn{h,T) 




det 


ds 


dxidxj 


(l + 0(l/n)) 


(9-42) 


where (xi,..., Xn) is the solution of the saddlepoint equation which minimizes 3?iS. 

Remark 9.1 The saddle point equation may have more than one minimal solution (ir). 

- in particular if ^ G M, there are two solutions, complex conjugate of each other. 

- in the multicut case, there are many saddlepoints with near-minimal action. 

In all cases, one needs to sum over all the saddle points. Let us call {x}/, the collection of saddle points. 
We have: 


Z„ 


E 


Cl 






e T 


(l + 0(l/n)) 


(9-43) 


Each saddle point {x}i corresponds to a particular minimal n-dimensional integration path in C"",noted Tj, 
and the coefficients Cj € h are such that: 

^ C/E/ (9-44) 

/ 

10 Solution of the saddlepoint equation 

We recall the saddle point equation: 


We introduce the function: 


Vi = l,...n, = 2- ^ _ 

rp . _ . 

vO 2 ^ 7 


T _ ^ . 1 

u;{x) 


n ^' X — Xn 

J=1 J 


(10-45) 


(10-46) 


in the large N limit, ui{x) is expected to tend toward the resolvent, at least in the case there is only 
one minimal saddle point. Indeed, the Xj’s are the position of the eigenvalues minimizing the action, 
i.e. the most probable positions of eigenvalues of M, and thus 110-461 should be close to 


10.1 Algebraic method 

Compute ijj‘^{x) -f —uj'{x)^ you hnd: 


a;^(x) -I- —u}'(x) 
n 


'X'. ^ __ ,__ 

—YY - - / - 

^ ^ ix — Xi){x — Xn) rfi ^ (x — 

1=1 7 = 1 ^ 1=1 '' '■> 


7^2 ^^ 


^The validity of the saddle point approximation is not proven rigorously for large number of variables. But here, 
we have many evidences that we can trust the results it gives. The asymptotics we are going to find have been proven 
rigorously by other methods. Basicaly, it is expected to work because the number of variables n is small compared to 
the large parameter in the action. 















n 


1 


(10 - 47) 



T y VL{x) - {VL{x) - 

n ^ X — Xi 

i=l * 

h . , ^ T ^V\x)-V'i:xi) 


{V'{x) - -z)uj{x) -^ 

X — ^ n “ 


i=l 




+ h 


^(0 

x-^ 


i.e. we get the equation: 

Uj^(x) + -Uj'(x) = v'(x)uj(x) - P(x) - (10-48) 

n X — ^ 

where P{x) := ^ ^ ^ polynomial in x of degree at most deg 17 — 2. 

In the large N limit, if we assume'^ that we can drop the l/NW'{x) term, we get an algebraic 
equation, which is in this case an hyperelliptical curve. In particular at h = 0 and T = 1: 

oj{x) = ^ ^17'(a;) — \/V'‘^{x) — 4P(a;)j (10-49) 

The properties of this algebraic equation have been studied by many authors, and the T and h 
derivatives, as well as other derivatives were computed in various works. Here, we briefly sketch the 
method. See [2nil2HllZI] for more details. 


10.2 Linear saddle point equation 

In the large N limit, both the average density of eigenvalues, and the density of x tend towards a 
continuous compact support density p(x). In that limit, the resolvent is given by: 

/X rj. f p{x')dx 

^{x) = T / - - 

_ J supp p 

is possible to do the calculation without droping the 1/iV term. One gets a Ricati equation, which is equivalent 
to a Schroedinger equation. If one is interested in a large N limit for the resolvent, the asymptotic analysis of that 
Schroedinger equation (Stokes phenomenon) gives, to leading order, the same thing as when one drops the 1/7V term. 
If one whishes to go beyond leading order, many subtleties occur. 


t 

(10-50) 



















i.e. 


'ix e supp p, p{x) = 


2iTiT 


[ijj{x + iO) — uj{x — iO)) 


and the saddle point equation 110-451 becomes a linear functional equation: 

Vx e supp p, + *0) + ~ *0) 


(10-51) 


(10-52) 


The advantage of that equation, is that it is linear in uj, and thus in p. The nonlinearity is hidden 
in supp p. 

10.2.1 Example: One cut 

If the support of p is a single interval: 


supp p = [a, 6] , a < b 


then, look for a solution of the form: 


= \ {Vhix) - Mh{x)^J{x-a){x-h)) 


(10-53) 


(10-54) 


The saddle point equation 110-521 implies that Mh{x + iQ) = M/i(x —iO), i.e. has no discontinuities, 
and because of its large x behaviour, as well as its behaviours near it must be a rational function 
of X, with a simple pole at x = ^. Mh, a and b are entirely determined by their behaviours near 
poles, i.e.: 

T 

(10-55) 


T 

uj{x) ~ — 

x^oo X 


ui{x) ~ regular — Mh{x) 


h 


Thus, one may write: 


a;(x) = - [ V'{x) — M(x)\/(x — a)(x — b) - 

2 \ x-^ 


X - ^ 

-^(x — a)(x — b) 

~ VK-«)(?-¥ 


(10-56) 


(10-57) 


where M(x) is now a polynomial (which still depends on h and T and the other parameters), it is 
such that: 


M(x) = Pol 


V'{x) 


-^(x — a)(x — b) 


(10-58) 


The density is thus: 


p{x) = l^Mh{x)^J{x -a){b- x) , supp p = [a, b] 


(10-59) 



















10.2.2 Multi-cut solution 

Let us assume that the support of p is made of s separated intervals: 

supp p = 


(10-60) 


then, for any sequence of integers ni, n 2 ,..., such that n, = n, it is possible to find a solution 
for the saddle point equation. That solution obeys 110-521 as well as the conditions: 

r^i 


Vz = l,...,s , 


rii 


' O'! 


p{x)dx = Tj^ 


(10-61) 


The solution of the saddle point equation can be described as follows: 
let the polynomial a{x) be dehned as: 


cr(x) := ]^(x - ai){x - h) 


(10-62) 


i=l 


The solution of the saddle point equation I10-52L is of the form: 

= {Vh{x) - Mh{x)^/^i^ 


(10-63) 

where Mh{x) is a rational function of x, with a simple pole at x = Mh-, and cr(x) are entirely 
determined by their behaviours near poles, i.e.: 

ui{x) 


ui{x) ^ regular 


T 

x—^oo X 

Mdx) 


h 


< x-i 


(10-64) 

(10-65) 


and by the conditions that: 


Mh{x)\/a{x)dx = 2inT— 


’ at 


Vi = l,...,s , 


( 10 - 66 ) 















10.3 Algebraic geometry: hyperelliptical curves 

Consider the curve given by: 

^ {yl^{x) - Mh{x)\/{x -a){x- b)^ (10-67) 

It has two sheets, i.e. for each x, there are two values of uj{x), depending on the choice of sign of the 
square-root. 

- In the physical sheet (choice +^), it behaves near oo like uj{x) ~ T/x 

- In the second sheet (choice —h behaves near oo like uj{x) 

Since a;(x) is a complex valued, analytical function of a cmplex variable x, the curve can be 
thought of as the embedding of a Riemann surface into C x C. 

I.e. we have a Riemann surface S, with two (monovalued) functions defined on it: p ^ S, —^ 
x{p) G C, and p & S, ^ uj{p) G C. For each x, there are two p E £ such that x(p) = x, and this is 
why there are two values of a;(x). 

Each of the two sheets is homeomorphic to the complex plane, cut along the segments [a*, bi], and 
the two sheets are glued together along the cuts. The complex plane, plus its point at infinity, is 
the Riemann sphere. Thus, our curve T, is obtained by taking two Riemann spheres, glued together 
along s circles. It is a genus s — 1 surface. 









































10.4 Genus zero case (one cut) 

If the curve as genus zero, it is homeomorphic to the Riemann sphere S = C. One can always choose 
a rational parametrization: 

+7(P+1/p) , = (10-68) 

\/ {x — a){x — h) = 7 (p — 1/p) (10-69) 

so that cn is a rational function of p. 

That representation maps the physical sheet onto the exterior of the unit circle, and the second 
sheet onto the interior of the unit circle. The unit circle is the image of the two sides of the cut 
[a, 6], and the branchpoints [a, 6] are maped to —1 and -|-1. Changing the sign of the square root is 
equivalent to changing p — 1/p. 

The branch points are of course the solutions of dx/dp = 0, i.e. dxij)) = 0: 

dx{p) = 7 ^1- dp , dx{p) = 0 p = ±1 <->• x{p) = a,b (10-70) 

There are two points at cx), p = cx) in the physical sheet, and p = 0 in the second sheet. 




Since the resolvent a;(p) is a rational function of p, it is then entirely determined by its behaviour 
near its poles, the poles are atp = oo, p = 0,p = P( and p = p^ (the two points of S such 
that x{p) = such that p^ is in the physical sheet, and p^ is in the second sheet): The boundary 
conditions: 


uj{p) 

uj{p) 

Uj{p) 

Uj{p) 


T 

_ 

p^oo x{p) 

V'{X{P)) 

h 


T 

x{p) 


x{p) 
~ regular 


h 

x{p) 


p^Pi 


< 


(10-71) 






10.4.1 T derivative 


Now, let us compute duj{p)/dT at x{p) fixed. Eq. 110-711 becomes: 

dijj{p) 1 

dT p^oo x{p) 
duj{p) 1 


dT x{p) 

du{p) 

—— ~ regular 
Ol 

duj{p) 

—— ~ regular 
ol 


(10-72) 


Moreover, we know that uj{x) has a square-root behaviour near a and b, in {x — a){x — b), and a 
and b depend on T, thus duj/dT may behave in ((a: — a){x — b))~^/‘^ near a and 6, i.e. doo/dT may 
have simple poles at p = ±1. 

Finaly, duj{p)/dT, has simple poles at p = 1 and p = —1, and vanishes at p = 0 and p = oo, the 
only possibility is: 


duj{p) 


dT 


p 


x{p) - 1 ) 

which is better written in terms of differential forms: 


1 dp 
p dx 


(10-73) 


du}{p) 

dT 


dx{p) 

x{p) 


dp 

p 


d Inp 


(10-74) 


the RHS is independent of the potential, it is universal. 
With the notation 18-341 we have: 


fl{p)dx{p) 


dp 

p 


A(p) = 


(10-75) 


10.4.2 h derivative 

The h derivative is computed in a very similar way. 

duj{p) 


oh P^oo 
doj{p) 1 

dh p^o x{p) 
duj{p) 1 


dh p^Pe x(p) — ^ 
duj{p) 

—— ~ regular 
dh p^p? 


(10-76) 




























implies that duj/dh can have poles at p = 
possibility is: 

du}{p) 


±1 and at p = and vanishes at p 
-PPi 


dh 


x{p) 


lip - Pi){.p‘^ - 1) 


i.e. 


duj{p) 

dh 


dx{p) 

x{p) 


which again is nniversal. 

With the notation 18-351 we have: 


dp dp 
P P-Pi 


din 


p 


P-Pi 


0. The only 


(10-77) 

(10-78) 


H{p,p^)dx{p) 


dp dp 
P P —- 

^ Pi 


H{Pi) = In 


Pi 


Pi-Pi 



(10-79) 


10.5 Higher genus 


For general genns, the cnrve can be parametrized by d-fnnctions. Like rational fnnctions for genns 
0, 6*-fnnctions are the bnilding blocks of fnnctions dehned on a compact Riemann surface, and any 
such function is entirely determined by its behaviour near its poles, as well as by its integrals around 
irreducible cycles. All the previous paragraph can be extended to that case. 

Let cx)+ and cx)_ be the points at inhnity, i.e. the two poles of x{p), with cx)+ in the physical 
sheet and oo- in the second sheet. Let p = Pi and p = pg be the two points of £ such that x{p) = 
and with p^ in the physical sheet, and p^ in the second sheet. 

The differential form ui{p)dx{p) is entirely determined by: 


a;(p)da;(p) ~ 

p^oo+ x[p) 

Uj{p)dx{p) ~ dV{x{p)) - 

x[p) x[p) 

uj{j))dx{j)) ~ 

p^Pi a;(p) - 4 
ui{j))dx{p) ~ regular 

p^Pi 




u{p)dx{p) = ^ 

n N 


Res uj{j))dx{p) = —T 


00 + 


Res u}{p)dx{p) = T + h 

oo_ 

Res ui{p)dx{p) = —h 

Pi 

Res ui{j))dx{p) = 0 

Pi 


(10-80) 


Since doj/dT^ h can diverge at most like {x — ai) near a branch point a*, and dx{p) has a zero at 
Oj, the differential form dudx/dT, h has no pole at the branch points. 


10.6 Introduction to algebraic geometry 

We introduce some basic concepts of algebraic geometry. We refer the reader to |221 for instance. 


















Theorem 10.1 Given two points qi and q 2 on the Riemann surface S, there exists a unique differ¬ 
ential form dSq^^q^ip), with only two simple poles, one at p = qi with residue +1 and one at p = q 2 
with residue —1, and which is normalized on the Ai cycles, i.e. 


' Res dSq^ q^ij)) = +1 
Res dSq^^q^ip) = -1 

\ p^q2 


\ 


dSq,,q2{p) = 0 


Ai 


dS is called an “abelian differential of the third kind”. 


(10-81) 


Starting from the behaviours near poles and irreducible cycles ll()-8()[ we easily hnd: 

duj{p)dx{p) 


H{p,Pi)dx{p) = 


kl{p)dx{p) = 
dui{p)dx{p) 


dT 


= -dS, 


00_p,00 


-ip) 


dh 


x{p) 


x{p) 

= -dSp^^oo.ip) = dSp^^oo+{p) -d\n{x{p) - x{p^)) 


(10-82) 


(10-83) 


Theorem 10.2 On an algebraic curve of genus g, there exist exactly g linearly independent “holo- 
morphic differential forms” (i.e. with no poles), dui{p), i = 1,... ,g. They can be chosen normalized 
as: 

® dujip) = 5ij (10-84) 

JAi 

For hyperelliptical surfaces, it is easy to see that if L{x) is a polynomial of degree at most g — l = s —2, 


the differential form 




=dx is regular at oo, at the branch points, and thus has no poles. 


y/Yli=iO-o.i)0-bi) 

And there are g linearly independent polynomials of degree at most g — 1. The irreducible cycles Ai 
is a contour surrounding [oj, bi] in the positive direction. 


Definition 10.1 The matrix of periods is defined by: 

Tij := ® duj{p) (10-85) 

JBi 

where the irreducible cycles Bi are chosen canonicaly conjugated to the Ai, i.e. Ai fl Bj = 5ij. In our 
hyperelliptical case, we choose Bi as a contour crossing [ai,bi] and [as,&s]. 

The matrix of periods is symmetric Tij = Tji, and its imaginary part is positive ‘ssTij >0. It 
encodes the complex structure of the curve. 


The holomrphic forms naturaly dehne an embedding of the curve into O: 












Definition 10.2 Given a base point go ^ we define the Abel map: 


u{p) = {ui{p),...,Ug{p)) , Ui{p):= dufip) 


( 10 - 86 ) 


where the integration path is ehosen so that it does not eross any Ai or Bi. 

Definition 10.3 Given a symmetrie matrix r of dimension g, such that ^Tij > 0, we define the 
9-function, from C® —C by: 

e{u,T) = (10-87) 

in&a 

It is an even entire function. For any rh G it satisfies: 


e{u+ Tn)= 9(u) , 9{u + rrh) = 9{u) 


( 10 - 88 ) 


Definition 10.4 The theta function vanishes on a codimension 1 submanifold ofC^, in particular, 
it vanishes at the odd half periods: 


z = +Tm 2 ^ , m 2 G IF , (m\mi) G 2Z -|- 1 —> 9{z) = 0 

For a given such odd half-period, we define the characteristic z 9-function: 


(10-89) 


9fiu) := e*™2“+0(M + f) 


(10-90) 


so that: 


9fiu + m)= 9fiu) , 9z{u + rm) = 


(10-91) 


9S) = 0 


(10-92) 


Definition 10.5 Given two points p, q in S, as well as a basepoint po E S and an odd half period z, 
we define the prime form E{p, q): 


E{p,q) := 


Oz{u{p) - u{q)) 


(10-93) 


where dhfip) is the holomorphic form: 


d9Au) 

dhfip) ;= ^ ^ dui{p) 


^ dui - ^ 

i=l * “=0 


(10-94) 



Theorem 10.3 The abelian differentials can be written: 


With these dehnitions, we have: 

B^{u{p) - u{(X)-)) a^(p) - ^( 00 +)) 

^ 9g{u{p) — u{oO+)) ’ ^ p^oo+ 6 '£(m(cx0+) — m(cX)_)) 

_ -M(oo-))gj-(M((X)+) -M(pg)) _ 9g{u{oo+) - u{oo_)) dhs{p^) 

9^{u{p^)-u{p^))9s{u{oo+)-u{ooJ)) 9^{u{p^)-u{oo+)Y dx{pf) 


(10-95) 


(10-96) 

(10-97) 


11 Asymptotics of orthogonal polynomials 

11.1 One-cut case 

In the one-cut case, (i.e. genus zero algebraic curve), and if 1/ is a real potential, there is only one 
dominant saddle point if ^ ^ [(^,b], and two conjugated dominant saddle points if a; G [a,b]. More 
generaly, there is a saddle point corresponding to each determination of p^ such that x{p^) = i.e. 
p^ and pg = 1/p^. The dominant saddle point is the one such that 3fJ(I4ff(p^) — ^(0) is minimal. 
The two cols have a contribution of the same order if: 


3?Wfr(p^) = 3?Wfr(p^) 


(11-98) 


i.e. if C, is such that: 




rPe 


W{x)dx = 0 




(11-99) 


If the potential is real, it is easy to see that the set of points which satisfy iri-991 is [a, 6 ], in general, 
it is a curve in the complex plane, going from a to b, we call it the cut [a, b] (similar curves were 
studied in EH). 

Then we have: 


• For X ^ [a, b], we write ^ ^ + 7 (p^ 1/p^), 7 = 

P»(0 ~ fmT) (A(Ps))”“" +0(1/JV)) 

i.e. _ 


PniO 


7 




{'yP^T~^ e-TT«ff(p?)-^(«))(l + 0(l/iV)) 


( 11 - 100 ) 

( 11 - 101 ) 













For X G [a, b], i.e. p is on the unit circle p = e*'^, ^ ^ + 27 cos 0: 


+ + 0{1/N)) (11-102) 

- /n ■ ! p{x)dx-{n-N + h(l){^) + a]{l + 0{l/N)) (11-103) 

A/2sin0(4) \ Ja ^ J 


i.e. we have an oscillatory behaviour 


P(x) 



11.2 Multi-cut case 


In the multicut case, in addition to having saddle-points corresponding to both determinantions of 
we have a saddle point for each hlling fraction conhguration ni,... ,ns with = xi. We 

write: 


(11-104) 


The saddle point corresponding to Filing fractions which differ by a few units, contribute to the 
same order, and thus cannot be neglected. One has to consider the sommation over filling fractions 

HO]. 

Thus, one has to consider the action of a saddle point as a function of the filling fractions. We 
leave as an exercise for the reader to prove that the derivatives of F are given by: 


dF 

dei 



W{x)dx 


(11-105) 





















and: 


d^F 

deidT 

d‘^F 

deidh 


= —2in{ui{oo+) — Ui{oo-)) 
= -2i7r{ui{p^) -Mi(oo+)) 


d‘^F 

——— = -2niTij 
oeiOej 


The last relation implies that ^F is a convex function of e, thus it has a unique minimum: 

dF 


F , 3? 


dei 


We write: 


__^ ^ 

2171 dei 


= 0 


, 0 e 


(11-106) 

(11-107) 

(11-108) 

(11-109) 

( 11 - 110 ) 


We thus have the Taylor expansion: 


h 

F{T,h,^ ~ F{l,0,F)-277rC\e-F) + {T-l)— + -{V,^{p^)-V{O) 

frp _ I\2a2p l2 

+ -(T- l)hlnA{p,) - - InHip,) 

—2i7r{e — e")V(T— e*) — 2i7r(T — 1)(T— e")*('u(cx)+) — u{oo-)) 
—2i7ih{e — Fy{u{p^) — u{oo+)) -|- ... 


( 11 - 111 ) 


Thus: 

Z 


/ 

^ ^_7V2ir(i,o,?‘)g^"(-T-hff-|(V;ft(p)-V(0)-^l^f^+(T-l)hlnA(p) + ^lnr/(p)j 

P=P^^P^ 

^ ^ ^i7v(n—N^yT(n—N€^)^2i7vNB(‘^—^^) 

n 

^iTTN{T—l)(n—N^Y{u{oQ^)—u{oQ—))^i'KNh(fi—N^)^{u(j})—u{oop^)) 


(11 - 112 ) 


In that last sum, because of convexity, only values of n which don’t differ from NF form more than 
a few units, contribute substantialy. Therefore, up to a non perturbative error (exponentialy small 
with iV), one can extend the sum over the nds to the whole Z®, and recognize a ^-function (see 

ITTHTfll : 

p=Pi,p^ 












——227rA^^(T—l)^*('u(c50_|_)—'u(oo_))g —227rA^^ h^''^ {u{p)—u{oo^)) 

6{N{( — re") + N{T — l)(-u(cxo+) — u{oo_)) + Nh(u{p) — u^oo^)), r) 

(11 - 113) 

with T — 1 = and h = 0 or h = 1/N, we get the asymptotics: 

x{p)=^ 

6{N{(- re*) + (n - N){u{oo+) - u{oo-)) + {u{p) - u{oo+)), r) 

6{N{( — T^) + {n — N){u{oo^) — -u(cx)_)), r) 

(11 - 114) 

Again, depending on we have to choose the determination of p^ which has the minimnm energy. 
If we are on a cnt, i.e. if condition 111-991 holds, both determinations contribnte. To snmmarize, 
ontside the cnts, the snm 111-1141 rednces to only one term, and along the cnts, the snm 111-1141 
contains two terms. 

12 Conclusion 

We have shown how the asymptotics of orthogonal polynomials (a notion related to integrability) is 
deeply related to algebraic geometry. This calcnlation can easily be extended to many generalizations, 
for mnlti-matrix models da cni El CHI, non-hermitean matrices {f3 = 1,4) [20], rational potentials 

0 , ••• 
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